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Abstract 

We compute the 2ra-point coupling constants in the high-temperature phase of the 
2d Ising model by using transfer- matrix techniques. This provides the first few terms 
of the expansion of the effective potential (Helmholtz free energy) and of the equation 
of state in terms of the renormalized magnetization. By means of a suitable parametric 
representation, we determine an analytic extension of these expansions, providing the 
equation of state in the whole critical region in the t, h plane. 



1 Introduction 



Despite its apparent simplicity and the fact that since more than fifty years the exact ex- 
pression of the free energy along the h = axis is exactly known, the two-dimensional Ising 
model still provides many interesting issues. Since the original Onsager solution Q, several 
other exact results have been obtained for this model. In particular, closed expressions for 
the two-point correlation function at h = and an S-matrix solution on the t = 0, h ^ 
axis [0] exist (for a review, see JJ). However, little is known for generic values of t and h and 
in particular, there is no exact expression for the free energy and for the critical equation of 
state in the whole (t, h) plane. 

In order to determine the Helmholtz free energy (also called effective potential) and the 
equation of state, we begin by computing the first terms of their expansion in powers of the 
magnetization in the high-temperature (HT) phase. The coefficients of this expansion are 
directly related to the n-point zero-momentum renormalized couplings g n that are also inter- 
esting in themselves, since they summarize relevant (zero-momentum) physical properties of 
the quantum field theory that describes the Ising model in the vicinity of the critical point. 
For this purpose we combine transfer matrix (TM) techniques, conformal field theory (CFT) 
methods, and general renormalization-group (RG) properties of critical systems. Recent 
works |3],[5| have shown that a combination of these approaches can lead to very accurate 
estimates of universal ratios in the 2d Ising model. In this respect this paper is the natural 
continuation of || in which the first nontrivial coupling was determined with the same 
techniques used here. 

Starting from the expansion of the free energy in powers of the magnetization in the HT 
phase, we determine approximate representations of the equation of state that are valid in 
the critical regime in the whole (t, h) plane. This requires an analytic continuation in the 
complex t-plane H0], extending the expansion valid for t > to the low-temperature phase 
t < 0. For this purpose, we use parametric representations |^-[T0[|, which implement in a 



rather simple way the known analytic properties of the equation of state (Griffiths' analyt- 
icity). We construct a systematic approximation scheme based on polynomial parametric 



representations J7| and on a global stationarity condition ||11|| . This approach was success- 
fully applied to the three-dimensional Ising model leading to an accurate determination of 
the critical equation of state and of the universal amplitude ratios that can be obtained from 
it 0,0]. 



The paper is organized as follows: In Sec. 2 we set up the formalism and define the 
quantities that parametrize the free energy for small values of the magnetization. These 
coefficients are evaluated numerically in Sec. 3 by using TM techniques. In Sec. 4 we use a 
parametric representation to analytically extend this expansion to the whole critical region, 
obtaining approximate, but still quite precise, expressions of the critical equation of state in 
the whole (t, h) plane. Finally in Sec. 5 we draw our conclusions. In App. A we report our 
notations for the universal amplitude ratios, in App. B we present a detailed discussion of 
the finite-size behaviour of the free energy and of its derivatives, and in App. C we explicitly 
compute eh, one of the coefficients appearing in the expansion of the nonlinear scaling field 
associated with the magnetic field. 
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2 Small-field expansion of the effective potential in the 
high-temperature phase 

In the theory of critical phenomena continuous phase transitions can be classified into uni- 
versality classes determined only by a few basic properties characterizing the system, such as 
space dimensionality, range of interaction, number of components and symmetry of the order 
parameter. RG theory predicts that, within a given universality class, the critical exponents 
and the scaling functions are the same for all systems. Here we consider the two-dimensional 
Ising universality class, which is characterized by a real order parameter and effective short- 
range interactions. A representative of this universality class is the standard square-lattice 
Ising model defined by the partition function 

2 = ^ e 13 ^^^ anam+h J2n a "^ (2.1) 

<Xi=±l 

where the field variable a n takes the values {±1}, n = (n ,ni) labels the sites of a square 
lattice of size L Q x L\ y and (n,m) denotes a lattice link connecting two nearest-neighbour 
sites. In our calculations we treat asymmetrically the two directions. We denote by no the 
"time" coordinate and by ri\ the "space" one. We indicate by iV = L L\ the number of sites 
of the lattice and define the reduced temperature 

tsh w-< < 2 - 2 > 

where (3 C = log (y/2 + l)/2 « 0.4406868 ... is the critical point. 
We introduce the Gibbs free-energy density 

F(t,h) = ^log(Z(t,h)), (2.3) 

and the related Helmholtz free-energy density — in the field-theoretical framework usually 
called effective potential — 

F(t,M) = Mh-F(t,h), (2.4) 

where M is the magnetization per site. For t > 0, F(t, h) and F{t, M) have regular expan- 
sions in even powers of h and M respectively. Explicitly 

F(t,h)-F(t,0) = J2 7^y,X2n(t)h 2n , (2.5) 

n=l l /n J- 

F(t, M) - F(t, 0) = T.j^xZ'm^, (2.6) 

where X2n{t) is the zero-momentum n-point function and Xzn {t) * s ^ s one-particle irreducible 
counterpart. Note that Xzit) = x{t>)i where x{t) is the magnetic susceptibility. 
Expansion (|2.6|) can also be written in the equivalent forms 

r(t,M)-T(t,0) = G ^(l^ + T-^-g 2nV A (2.7) 



l& 2+ r 4+ S(4 r - 2 ")- (2 - 8) 
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where 



^ = ^^!M 2 , z 2 = -^#M 2 , (2.9) 

and ^2nd is the second- moment correlation length 

6nd = ^E^(^(0)^)) • (2-10) 

The coefficients g 2 „ correspond to the 2n-point renormalized coupling constants and r 2n = 
92n/g2~ 1 - The couplings #2n and the related r 2n can be computed in terms of \2n and ^2nd- 
The four-point coupling is given by 

9*(t) = (2.H) 

A S2nd 

An explicit expression of T2n in terms of \2n is given at the beginning of Sec. |||. 

The advantage of the expansions (|2.7[) and (|2.8|) is the fact that the coefficients have a 
finite limit — in the following we indicate it by the same symbol — for t — > 0. Notice that for 
t — > 0, z ~ Mt~P and X 2 /x<i ~ t 2u i so that Eq. ( |2.8| ) defines the HT expansion of the scaling 
part of the free energy F(t, M) ~ t 2u F acid (Mt-P). 



3 Calculation of rm with transfer- matrix techniques 

In this section we wish to report the calculation of the first few coefficients r2„ at the critical 
point t — 0. It is easy to express them in terms of the critical amplitudes of the 2n-point 
functions X2n{t)- Defining 

G„ + = lim Xn (t)t 15 "/8-2 ; (31) 
it is possible to show by direct calculation that 



r 6 = 10 (3-2) 

s~i+/~i+ ri+ ( y-f+\2 

r 8 = 280 - 56—^-+ , (3.3) 

n+n-v r<+(r i +\ 2 



r 10 = 15400 - 4620^ + 120- 



y 4 



1 io«(^6 h ) 2 ( ( ^2 h ) 2 Clo(C2y ,,, 
+ 126 ^^? (^"' (3 ' 4) 

r 12 = 1401400 - 560560^^ + 17160^^- + 36036 ( ° 6 j ( ° 2 } 



(ct) 2 (ct) 3 {cif 

~ 220 ^cfy~ ~ 792 (c+)* + "IcfF" (3 ' 5) 

The aim of this section is to estimate the constants C^n, by using TM techniques and the 
exact knowledge of several terms of the small-t expansion of X2n{t)- Using the equations 
reported above, we will then obtain our estimates of r 2ri . 
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3.1 The transfer-matrix approach 

We obtained our numerical estimates of the 2n-point functions X2n(t) following a two-step 
procedure. First, by using TM techniques, we obtained estimates of these quantities on 
lattices of finite transverse size Li, with L\ < 24. Second, we extrapolated these results to 
the thermodynamic limit by using the exact knowledge of their Finite Size Scaling (FSS in 
the following) behaviour. To perform this second part of the analysis we have elaborated a 
new extrapolation scheme which is rather interesting in itself. We shall discuss it in detail 
in the last part of this section. 

Let us now see both these steps in more detail. 

3.1.1 The TM computation 

This part of the procedure was already discussed in Ref . [[J . We report here the main features 
of the algorithm for completeness and refer to || for a more detailed discussion. 

The main idea is to use TM techniques to extract the /i-dependence of the magnetization 
at fixed t. To this end, we computed numerically the magnetization M(t, h, L\) on lattices 
of size oo x Ii, for L\ < 24. In this work we use again the data obtained in Ref. ||. In 
addition, we generated new data for (3 = 0.335 and 0.345 and added L\ = 24 results for 
(3 = 0.35 and 0.355. All computations were performed with double-precision floating-point 
arithmetic on commercial workstations. Typically we obtained M with 15 correct digits. 

Then, for each given value of L\ and (3, we determined the coefficients of the series 

h = btM + b 3 M 3 + b 5 M 5 . . . (3.6) 

which we truncated at order M 15 . Using the numerical results of M for 8 values of h, we 
computed the coefficients bi, 63, 615. The optimal choice for these eight values of h is 
the one for which the errors due to the truncation of the series and those due to numerical 
rounding are of the same magnitude. This optimal range changes as a function of j3 and L\. 
For instance, just to give an idea of the magnitude of the magnetic fields which we studied, 
for (3 = 0.36 and L\ = 24 the best choice is h = 0.002 j with j — 1 . . . 8 . The accuracy of 
the bi obtained in this way is decreasing with increasing order. For example, for (3 = 0.37, 
we obtain b\ with 14 significant digits and, e.g., bn with only 3 significant digits. 
In principle we could have also used the expansion 

M = a 1 h + a 3 h 3 + a 5 h 5 . . . . (3.7) 

However, in practice, the use of (|3.7| ) leads to results with significantly larger truncation 
errors, since in contrast to what happens for bi, the sign of the coefficients alternates. 

Then, for each (3 and Li, we constructed the X2n functions out of the bi constants. These 
are the inputs of the FSS analysis discussed in the next section. The relations between the 
bi and the X2n can be easily obtained by a straightforward calculation. We report them here 
for completeness: 



Xl2 



Xio 




(3.10) 



(3.11) 



(3.12) 



This last step is the only difference in the present calculation with respect to the analogous 
one reported in || where we directly studied the thermodynamic limit of the h\ coefficients. 
The reason is that as n increases the \ 2n show a much better and smoother FSS behaviour 
than the &j. 

3.1.2 The thermodynamic limit 

In the second step of the analysis we extrapolate our results to the thermodynamic limit. We 
know on general grounds that, since L\ ^> £, the convergence towards the thermodynamic 
limit of the X2n functions is exponentially fast, the decay rates being related to the spectrum 
of the theory. This result holds for any lattice model. However, in the case of the two 
dimensional Ising model, thanks to the fact that the model can be solved exactly also on 
finite lattices ]12] (or, in the language of S-matrix theory, thanks to the fact that the S- 
matrix of the model is very simple), much more informations can be obtained on the FSS 
properties of the free energy and of its derivatives. In particular, one can explicitly compute 
the functional form of the FSS behaviour of \2n (see App. [B]). It turns out to be 



X2n{L 1 ) -X2n(oo) = Lf"^ (g 2nA (l/L^e'™^ + g 2n , 2 {1/L X ) e" 2mLl + g 2n , 3 {1/L X ) e" 3mLl + . . .) 



where g 2n ,i(^/Li) are functions of 1/Li, which can be expanded in (positive) powers of 

By fitting our values of X2n{L\) with this law (by expanding the functions g2n,i{Li) and 
using as free parameters the coefficients of the first few terms of the expansion) it is possible 
to obtain very precise estimates for the thermodynamic limit X2n(oo). However, this is not 
the most efficient strategy, since much of the information contained in the X2n(Li) is lost in 
the determination of the coefficients of the Taylor expansion of the g2n,i{Li) functions. Thus, 
we have elaborated an alternative iterative procedure which is much simpler to perform than 
the above multiple fits and allows to reach higher precision (up to one additional accurate 
digit). The idea behind this Iterative Algorithm (IA, in the following) is first to absorb 
the pre-exponential factors of Eq. ( |3.13| ) in the masses by allowing them to depend on L\. 
Then, the resulting exponential corrections are eliminated by iteratively solving the system 
of equations 



(3.13) 



" 2) 
xHl{U - 1) 



c exp(— x(Li 
c exp(—x(Li 
c exp(— xL\) 



2)) + 

1)) + xt^iLi) 

+ xH: x \li) 



(3.14) 



with respect to xin C^i)? c an d x 
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The index i in xin denotes the order of the iteration, in particular the X2I are the input 
data of the algorithm obtained as discussed in the previous section. As i increases the L\ 
dependence of the x£n becomes smaller. After a certain number of steps the extrapolation 
becomes unstable since rounding errors accumulate. Typically, the final estimate was ob- 
tained by iterating 1, 2, 3 or 4 times, depending on (3 and on n. The residual dependence 
on Li is used to estimate the error in the best estimate of X2n(oo). The results obtained 
in this way are in perfect agreement with those obtained by directly fitting Eq. ( [3.13|) but, 



as mentioned above, turn out to be more precise. By using the IA we were able to obtain 
reliable estimates of \n for n < 12 up to (3 = 0.37. Our results are summarized in Tables [I] 
and H 

In order to verify if our method of estimating the errors with the IA method (quoted in 
Tables [I] and 0) is reliable we made the following test. We studied with the IA a sample of 
data at a very low value of (3 where rather large values of could be reached, so that the 
results for the X2n{Li) on our largest lattices were already good estimates of X2n(oo) with 
no need of further manipulations. Then, we performed our analysis using only the data at 
small L\. In this way we could explicitly check that our extrapolation scheme gives accurate 
results for the thermodynamic limit starting from data such that Li/£ ~ 7, which is what 
we reached at our largest value of /3, (3 — 0.37. Moreover in App. p.3| we report a further 
test of IA. We applied it to a test function of the type (|3.13|) , finding again support to the 
reliability of the method. 

In order to give a feeling of the performances of the IA for j3 = 0.37 — this is the value of 
13 that is nearer to the critical point among those we studied and hence the worst one from 
the point of view of FSS — we have reported in Table |3] the results of the first four iterations 
for %4 at /3 = 0.37. From these numbers we estimate the thermodynamic limit of %4 to 
be ^4(00) = —4.052238(1) to be compared with the result obtained by directly fitting Eq. 
( frig ) which is x 4 (oo) = -4.052242(10). 

This same IA was also used in our previous paper |J, the only difference being that in 
that case it was applied to extract the thermodynamic limit of the hi constants. Let us 
stress, as a final remark, that this algorithm is very general and can be used even when no 
exact information is available on the FSS behaviour of the quantity of interest except from 
the fact that it is dominated by an exponential decay. 

3.2 Small-£ expansion of Xn(t) 

As it is well known, the TM approach gives reliable results only for rather small values of 
(3. In order to perform the extrapolation (3 — > /3 C , it is thus mandatory to have a good 
control of the scaling corrections of Xn(t)- in this section we address this problem in detail. 
Our goal will be to obtain for each \in the exact form of the scaling function up to 0(t 4 ) 
and the spectrum of the possible scaling dimension (i.e. the exponents of the corresponding 
contributions in the scaling function) up to 0(t 5 )). 

The most important information which is needed to perform this analysis is the spectrum 
of the irrelevant operators of the theory. This problem was recently addressed numerically 



in Refs. ||i3|-4T5f where it was shown that, for rotationally invariant quantities like the free 
energy and its derivatives, the first correction due to the irrelevant operators in the square- 
lattice Ising model appears at order t A (hence no correction of order t 2 is present) and that 



6 



p 


4 


6 


0.20 


-3.2111149829(1) 


61.94504968(2) 


0.25 


-3.4721394468(1) 


74.5437488(1) 


0.28 


-3.622540033(1) 


82.255117(3) 


0.30 


-3.720514859(2) 


87.459547(5) 


0.31 


-3.768883189(5) 


90.081633(5) 


0.32 


-3.81687386(2) 


92.71755(1) 


0.33 


-3.86451567(5) 


95.36788(5) 


0.335 


-3.88821530(5) 


96.6986(1) 


0.34 


-3.91183945(10) 


98.0332(1) 


0.345 


-3.9353922(1) 


99.3718(1) 


0.35 


-3.9588780(5) 


100.7142(2) 


0.355 


-3.9823005(10) 


102.0610(2) 


0.36 


-4.0056653(10) 


103.4119(2) 


0.365 


-4.028975(1) 


104.7673(2) 


0.37 


-4.052238(1) 


106.127(1) 



Table 1: Thermodynamic-limit results for Xn(t) t 15n / 8 2 , for n = 4, 6. The quoted error bars 
are estimates of the systematic error of the extrapolation. 



P 


8 


10 


12 


0.20 


-2985 


5224(2) 


267883. 


:io 


-23211300. 


(1000.) 


0.25 


-4002 


0604(2) 


400060. 




-43255500. 


(1000.) 


0.28 


-4672 


309(2) 


494203. 




-59550000. 


(10000.) 


0.30 


-5144 


20(2) 


563455. 


10.) 


-72400000. 


(10000.) 


0.31 


-5387 


715(10) 


600090. 


;2o.) 


-79445000. 


(30000.) 


0.32 


-5636 


31(2) 


638110. 


;5o.) 


-86910000. 


(100000.) 


0.33 


-5890 


10(5) 


677600. 


100.) 


-94830000. 


(100000.) 


0.335 


-6018 


9(1) 


697800. 


100.) 


-98950000. 


(200000.) 


0.34 


-6149 


0(1) 


718300. 


100.) 


-103190000. 


(200000.) 


0.345 


-6280 


5(2) 


739000. 


;5oo.) 


-107560000. 


(200000.) 


0.35 


-6413 


1(3) 


760600. 


;2oo.) 


-112000000. 


(500000.) 


0.355 


-6547 


1(3) 


782400. 


;2oo.) 


-116600000. 


(500000.) 


0.36 


-6682 


3(5) 


804400. 


;3oo.) 


-121300000. 


(1000000.) 


0.365 


-6819 


3(3) 


827100. 


;5oo.) 


-126100000. 


(1000000.) 


0.37 


-6958 


1(5) 


850000. 


;5oo.) 


-131000000. 


(1000000.) 



Table 2: Thermodynamic-limit results for Xn(t) t 15n ^ 8 2 for n = 8, 10, 12. The quoted error 
bars are estimates of the systematic error of the extrapolation. 
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Li 


fP 




fP 




fP 


15 


-3.72609418989519 


-4.07455655 








16 


-3.80380147054199 


-4.06509498 








17 


-3.86341299292473 


-4.05978684 


-4.05300296 






18 


-3.90897624823899 


-4.05675321 


-4.05270713 






19 


-3.94370300243699 


-4.05498869 


-4.05253535 


-4.05229750 




20 


-3.97011104914228 


-4.05394438 


-4.05243012 


-4.05226371 




21 


-3.99015681870086 


-4.05331548 


-4.05236336 


-4.05224754 


-4.05223270 


22 


-4.00535064954957 


-4.05293008 


-4.05232014 


-4.05224075 


-4.05223585 


23 


-4.01685278847850 


-4.05268980 


-4.05229191 


-4.05223874 


-4.05223789 


24 


-4.02555119351443 


-4.05253744 


-4.05227342 


-4.05223837 


-4.05223829 



Table 3: Results of the IA analysis for F A = 11//2 at f3 = 0.37. 

the next irrelevant operator contributes only at order t G . This result has been shown for the 
susceptibility in zero field, but a standard application of RG ideas implies that it should also 
apply to the free energy as a function of t and This is also confirmed by the less precise 
numerical results of Ref. |4[] for the free energy on the critical isotherm and by the analytic 
study of the two-point function at large distances JKJ]. 

This result allows us to determine the scaling corrections to X2n(t) up to order t 4 by using 
standard RG techniques (see |5|, |T7| , [T8| ]). 

As a first step, let us write the free energy of the model in terms of nonlinear scaling 
fields |]| : 



\d/vt f . I i U 3 



F(t,h) = F b (t,h) + \u t \ d ^f smgu . , 




+N*logkl/..n 5 (^,{^n- (3.15) 

Here Fb(t,h) is a regular function of t and h 2 , Ut, Uh, {uj} are the nonlinear scaling fields 
associated respectively to the temperature, the magnetic field and the irrelevant operators, 
and y t , yt, {yj} are the corresponding dimensions. For the Ising model y t = 1 and y^ = 15/8. 
Notice the presence of the logarithmic term, that is related to a "resonance" between the 
thermal and the identity operator.^ Since all numerical data indicate that for t — > all zero- 
momentum correlation functions diverge as a power of t without logarithms — our results 
provide additional evidence for this cancellation — we shall assume in the following (as in 

lr This pattern agrees with an independent analysis, in the framework of CFT, based on the spectrum of 
the quasi-primary operators of the critical Ising model. This result was already anticipated in Ref. [^| and 
will be discussed in full detail in a forthcoming publication. 

2 In principle, other logarithmic terms may arise from additional resonances due to the fact that yj are 
integers or differ by integers from yt, and indeed they have been observed in a high-precision analysis of the 
asymptotic behaviour of the susceptibility |]l5| . They will not be considered here since these contributions 
either are subleading with respect to those we are interested in or have a form that is already included. 
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Ref. |17[]) that / s i ng does not depend on its first argument Uh/\u t \ yh ^ Vt . There is also some 



evidence that the leading contribution due to the irrelevant operators is absent. Indeed, 
for the susceptibility one would expect a correction of order t 4 log \ t\ which is not found in 
the high-precision study of the susceptibility reported in [15]. In the following we will be 
conservative and we will report results with and without corrections of order t 4 log Our 
final estimates assume however that such a term is absent. 

The scaling fields are analytic functions of t and h that respect the Z 2 parity of h. Let 
us write their Taylor expansion, keeping only those terms that are needed for our analysis 
(we use the notations of Ref. [jT^[ ): 

(3.16) 
(3.17) 



u h = h [1 + c h t + d h t 2 + e h h 2 + f h t 3 + 0(t\th 2 )}, 
u t = t + b t h 2 + c t t 2 + d t t 3 + e t th 2 + f t t 4 +0(t 5 ,t 2 h 2 ,h*). 



All these coefficients (except e/J have been determined exactly or numerically with very high 
precision, see &|I3,[18[]. We list them here for completeness: 



Ch 



y/2' 



— 1=, dt 



23ff 
16 

6 ' 



fh 
ft- 



191/3 C 3 
6V2' 



btPcV2, 



EqTT 

m 



2 : 

c 



where |0] 



E = 0.0403255003... 



(3.18) 

(3.19) 
(3.20) 

(3.21) 



is the coefficient of the contribution proportional to t log |t| in the susceptibility. 

The coefficient is the only term that was not reported in [[J. In App. y we show that 
it is possible to determine it by analyzing the scaling corrections to the free energy on the 
critical isotherm. Using the precise data of Ref. H, we obtain 



e h = -0.00727(15). 



(3.22) 



Using Eq. ( 3.15|) and the expansions of the scaling fields reported above, we can compute 
the leading terms in the asymptotic expansion of Xnif) f° r t — > 0. We obtain 



where 



Xn (t) = t 2 - 1 W8 an(t)+t 19/4-15n/8 6n(t) + 



a n (t) = C+il + a^t + azP^ + a^ + Oit 4 )), 
Ut) = Ct 2 ((o + (it + 0(t 2 )). 



The coefficients ati and Q are given by: 

7n - 16 



a i 



8V2 



(3.23) 



(3.24) 
(3.25) 



(3.26) 
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147n 2 - 1080n + 2176 
°* = Wg , (3-27) 

343n 3 - 5208n 2 + 26240n - 49152 

«3 = 7= , 3.28 

6144v/2 V ; 

E n(n-l)(15n-46)7r 
Co - ^ , (3.29) 

, E n(n - l)(7n - 38)(15n - 46)vr 

Ci = +n(n-2)(n-lK (3.30) 

Plugging these coefficients into Eq. ( |3.23| ) we obtain the exact form of the scaling function 
up to the contribution of the first irrelevant field, i.e. up to 0{t A ). The simplest way to use 
the exact knowledge of these terms of the scaling function to extract the amplitudes C + is 
to construct the quantity: 

V (t) f!W8-2 _ U / f \ f 15n/8-19/4 
M (n) (t) = XnWl WW (3.31) 

a n (t) 

which has the following expansion for t — > 0: 

M<">(t) =C+(l+ Pl t* + p l t i logt + p 2 t i - 75 + 0(t 5 )). (3.32) 

The terms proportional to t 4 and to t 4,75 are due to the first unknown coefficients of a n (t) and 
b n (t) respectively. They also take into account the possible presence of irrelevant operators 
contributing to order t 4 . The term proportional to t 4 logi is the only remnant of the / S i ng 
term in Eq. fl3.15| ) and it is due to the irrelevant operators in / S i ng . 



The constants C + are determined in sequence. We start by using C%, which is known 
to very high precision, C 2 + = 0.9625817323087721140443 . . . to estimate C 4 + with 

a best-fit analysis of Eq. Q3.32 ) with n — 4. The value of C£ obtained in this way is then 



used as input to construct the function 6 6 (£) which appears in Eq. fl3.31| ), thus allowing to 
estimate M^{t). At this point, by using again Eq. ( |3.32| ), we obtain the best-fit value for 
Cq and repeat the whole procedure for the next value of n. 

As a consequence, in the determination of (t) there are three sources of uncertainty: 

a] The uncertainty in the TM estimates of Xnif) which we use as input of our analysis. 

b] The uncertainty in the estimate of and Eq. 

c] The uncertainty in our best estimate of C^_ 2 . 

These uncertainties must be treated in different ways, [a] can be straightforwardly prop- 
agated to M^ n \t). On the contrary [b] and [c] appear as uncertainties of the whole fitting 
function. To deal with them we followed the simplest (and most conservative) strategy. Let 
us study as an example the uncertainty due to the error on e^. We constructed two sets of 
data, M+\t) and M_(t), obtained by using eh + Set and e^ — Set respectively, and then, for 
each set, we determined the best-fit values of the parameters. The difference between these 
results is a conservative estimate of the error induced by the uncertainty on e^. In order to 
give an idea of the size of these uncertainties, we have reported in Table |5] those induced by 
et in the best-fit estimates of Cq. 
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It turns out that the errors induced by the uncertainties on C+_ 2 and E are always 
negligible with respect to that due to and, what is more important, that this last one is in 
any case negligible with respect to the uncertainty in the best-fit estimate, i.e. with respect 
to the uncertainty due to the lack of knowledge of the higher-order correction terms. For 
this reason we shall neglect the errors of type [b] and [c] in our final results. 

3.3 The fitting procedure 

To analyze the data we followed a procedure similar to that presented in Ref. ||. For each 
set of data we performed five different fits of using Eq. (|3.32| ), 

fl] keeping and p\ as free parameters, and setting pi = p 2 = 0; 

f2] keeping and p\ as free parameters, and setting pi = p 2 = 0; 

f3] keeping C+, pi, and pi as free parameters, and setting p 2 = 0; 

f4] keeping C+, p±, and p 2 as free parameters, and setting pi — 0; 

f5] keeping C+, pi, pi, and p 2 as free parameters. 

In order to estimate the systematic errors involved in the determination of C+, we per- 
formed for all the fitting functions several independent fits, trying first to fit all the existing 
data (reported in Tables [I] and |2|) and then eliminating the data one by one, starting from 
the farthest from the critical point. Among the set of estimates of the critical amplitudes 
we selected only those fulfilling the following requirements: 

1] The reduced \ 2 °f t ne fit must be of order unity. In order to fix precisely a threshold, 
we required the fit to have a confidence level larger than 30%. 

2] For all the subleading terms included in the fitting function, the amplitude estimated 
from the fit must be larger than the corresponding error. 

3] If the fit contains k free parameters besides C+, then at least 2 k degrees of freedom 
must be used in the fit. This means that in fits of type fl and f2 at least four data 
must be used {k = 1: two parameters plus two degrees of freedom) in those of type 
f3 and f4 at least 7 data must be used {k = 2: three parameters plus four degrees 
of freedom). In f5 at least 12 data must be used [k = 3: four parameters plus eight 
degrees of freedom) 

This last constraint is a generalization to higher values of n of the one that we proposed 
in ||. It seems to encode very well the real statistical significance of the data. 

Finally, among all the estimates of the critical amplitude fulfilling these requirements 
we selected the smallest and the largest ones as lower and upper bounds. We consider the 
mean of these two values as our best prediction for C+. The values are reported in Table [|. 
The errors quoted in Table |] are half of the difference between the upper and lower bounds. 
They seem to give a reliable estimate of the uncertainty of our results. 

Let us briefly comment on these results: 
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Pi 7^0 


Pi = 


ct 


-4.379095(8) 


-4.379094(6) 




125.9330(11) 


125.9332(6) 


c i 


-9066.5(9) 


-9066.4(7) 




1216.33(80) x 10 3 


1216.34(60) x 10 3 


C\2 


-26260(150) x 10 4 


-26175(60) x 10 4 



Table 4: Best-fit estimates of the amplitudes ■ 

a] In Table ^ we have reported in two separate columns the best-fit results with and without 

the log-type contribution pit A logt. The reason of this choice is that there are strong 
numerical indications JHJ that p\ = in the Ising model. To keep our analysis as 
general as possible, we report also the result with pi ^ 0. It is interesting to note that 
the estimates of (n < 10) do not depend essentially on this choice, while the error 
decreases by a factor 1.5-2 if we assume pi = 0. This seems to indicate that the fitting 
procedure is stable and reliable. 

b] In the first line of Table || we report our best-fit estimate of C£, which we had already 

estimated in || with the same techniques used in the present paper. The only improve- 
ment with respect to [[J is that we now also know exactly the value of eh which was the 
dominant (of order t 3 75 ) unknown correction in ||. It is instructive to compare the 
two estimates. In || we obtained, keeping also into account the log-type corrections, 
Ct = -4.379101(9), while the present value is C 4 + = -4.379095(8). The two results 
are, as expected, fully compatible. The new one is only slightly more precise than the 
previous one. The fact that the enhancement in precision is so small is due to the fact 
that, by using the exact value of eh, we only improve the known small-t expansion of 
Xn (t) from 0(t 3 - 75 ) to 0{t A ). 

c] It is instructive to look in more detail to the results of the fits in one particular case. 

Let us consider, for instance, C§ . The results are reported in Table |], where, for each 
type of fit, we quote the two fits which correspond to the highest and lowest values of 
Cq (in quoting the uncertainties in the best-fit estimates we also report those induced 
by eh, to give an idea of their size). 

The most impressive feature of this set of fits is that with only the t A contribution 
besides Cg" (fit fl) one can fit all the data up to (3 = 0.31. 

Looking at Table [|, we obtain our best estimate 

C 6 + = 125.9330(11). (3.33) 

If we additionally assume that p\ = 0, i.e. we keep into account only the results of the 
fits of type fl and f4, then we find 

Q + = 125.9332(6). (3.34) 

The estimates ( |3.33| ) and ( |3.34D are those reported in Table [|. 
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d.o.f. 


fit type 


125.93274(12)(0) 


2 


fl 


125.93298(5)(0) 


9 


fl 


125.93301(14)(2) 


2 


f2 


125.93389(16)(7) 


7 


f2 


125.93223(30)(1) 


5 


f3 


125.93341(15)(2) 


8 


f3 


125.93290(14)(0) 


7 


f4 


125.93369(12)(1) 


9 


f4 


125.93253(20)(2) 


8 


f5 


125.93258(7)(2) 


10 


f5 



Table 5: Estimates of Cq . In the first column we report the best- fit results for the critical 
amplitude (the first error in parenthesis is that induced by the systematic errors of the input 
data while the second is due to the error on eh), in the second column the number of degrees 
of freedom (i.e. the number of data used in the fit minus the number of free parameters) 
and in the last column the type of fit. 



3.4 Summary of the results 

Using the estimates of obtained in the previous section, we can estimate r 2n by using 
Eqs. (pD-(pD- We finally obtain 

#4 = 14.697323(20), (3.35) 

r 6 = 3.67866(3)(2), (3.36) 

r 8 = 26.041(8) (3), (3.37) 

no = 284.5(1.4)(1.0), (3.38) 

Ha = 4200(320)(420). (3.39) 

For each r 2n we report two errors: the first one is due to the error on while the second 
one expresses the uncertainty due to the error on all C^, k < n. Note that the high 
precision reached in the estimates of C% n is partially lost in the estimates of r 2n because of 
the cancellations among the various terms in the sums. This effect is particularly important 
for no and ri 2 . 

Tables ^] and compare our results with the existing data from other approaches. Our 
estimates of r 2n perfectly agree with those obtained analyzing HT expansions. On the 
contrary, the Monte Carlo results for r 2n of Ref. E| are systematically larger. This could 



be an indication that the finite-size scaling curves obtained in Ref. |25[ are still affected by 
large scaling corrections. 



4 The critical equation of state 
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Method 


Ref. 


#4 


TM+CFT 


[this work] 


14.697323(20) 


FF 


PI 




14.6975(1) 


HT 


22 




14.694(2) 


HT 


23 




14.693(4) 


HT 


24 




14.700(17) 


MC 


21 




14.69(2) 


MC 


[25 




14.7(2) 


FT e-expansion 


26 




14.7(4) 


FT g-expansion 


27 




15.5(8) 


FT g-expansion 


p8 




15.4(3) 


c?-expansion 


29 




14.88(17) 



Table 6: Estimates of g±. We also report the existing results from the form- factor approach 
(FF), high-temperature expansions (HT), Monte Carlo simulations (MC), field theory (FT) 
based on the e-expansion and the fixed-dimension d = 2 g-expansion, and a method based 
on a dimensional expansion around d = (d-expansion) . 





TM+CFT 


HT 


ht m 


MC H 


FT e-exp. [30 


FT g-exp. 31] 




3.67866(5) 


3.678(2) 


3.679(8) 


3.93(12) 


3.69(4) 


3.68 


r 8 


26.041(11) 


26.0(2) 




28.0(1.6) 


26.4(1.0) 




t*io 


284.5(2.4) 


275(15) 










7*12 


4200(740) 













Table 7: Estimates ofrm- We also report the existing results from high-temperature expan- 
sions (HT), Monte Carlo simulations (MC), and field theory (FT) based on the e-expansion 
and the fixed-dimension d = 2 g-expansion. 



4.1 General features 

The basic result of the RG theory is that asymptotically close to the critical point the 



equation of state may be written in the scaling form [32 



h = ^^M\M\ s - 1 f(x), xoztm- 1 ", (4.1) 

where f(x) is a universal scaling function normalized in such a way that /(— 1) = and 
/(0) = 1. The value x = — 1 corresponds to the coexistence curve, and x = to the critical 
point t — 0. The function h(M,t), representing the external field in the critical equation of 
state, satisfies Griffiths' analyticity, i.e. it is regular at t — for M > fixed and at M = 
for t > fixed. This implies that f(x) is analytic at x — 0, and it has a regular expansion 
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for large- x of the form 



/(*) 



x 



n=0 



(4.2) 



As already mentioned in the introduction, many things are exactly known for the two- 
dimensional Ising model. However, there is no exact expression for the free energy and for 
the critical equation of state in the whole (t, h) plane. In Tabled we report a summary of the 
known results for the two-dimensional Ising model (there we consider only infinite-volume 
quantities). Many of them are known exactly, for the others we report their best estimate. 
The results that have not been derived in this paper have been taken from Refs. [B|,fJ3-£JB|. 

In the following we will determine the equation of state, starting from its expansion 
for small magnetization in the HT phase. It is therefore useful to introduce a different 
representation that is analytic for M — > 0. Using the results of Sees. |2| and |3.2| , in particular 
Eq. ([3.15 ) and the discussion following it, one may write the Helmholtz free energy as 



^ sing (t, M) 



at 2 V(z) + -t 2 \og\t\ 
2 



where 



z = b\M\t-P, 



a = -{c+y/c 



2 m+ 

4 ) 



b = 



z 2 Z A 

y + 4T + ^ 6 



(4.3) 



(4.4) 



CI 



(c+y 



1/2 



The constant A is related to the amplitudes of the specific heat for h — > defined in Eq. ( |A.1| ). 
Indeed the analyticity of the free energy for t = 0, h ^ implies A + = A~ = A, and thus 
Uq = A + / A~ = 1. The presence of the logarithmic term gives rise to logarithms in the 
expansions of V(z) for z — > oo. Indeed, the analyticity of ^ing(t, M) for t — 0, \M\ ^ 
implies, for large z, 



V(z) 



.16 



E CkZ 8k + ciog log Z, 



k=0 



The constant cw is easily expressed in terms of invariant amplitude ratios: 



AA 

ciog = — = ^Q + gt- 

a 



For the equation of state we have 



h 



dM 



where, using Eq. ( |2.8| 



Biz) 



abt 15 / 8 



z + ±z 3 



dV(z) 
dz 



,bt 15/8 B{z) 



r 2j 



5 (2J " 1)! 



,2i-l 



For large z, using (|4.5|) we obtain the expansion 



B(z) 



. 15 



(4.5) 

(4.6) 

(4.7) 
(4.8) 

(4.9) 



fc=0 
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Critical exponents and amplitude ratios 
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7/4 


V 




1 


7/ 




1 1 A 

1/4 


P 




1/8 


5 




15 


UJ 




2 


T T 

Uq = 


A + 1 A — 

A^ / A 


-1 
1 


T T 

U 2 = 




6( .byooozUl 


D+ — 


: A L<2 1 D 


u.oioooyoy 


K c = 


■ A-C 2 /B 2 


U.UUo4olo4 


J) 

K x = 


Q l = Li B° L /{dC c )° 


0. ( ( ozooUz 


„ ,2 

W = 


/-»— l\ r>2( t—\2\ 

C 2 /[B {J ) \ 


0.53152607 


Ut = 




3.16249504 


T J 


= f + If" 

J gap/ ./gap 
- A + 1 -f +\2 


2 


Q + = 


0.15902704 


Q = 


A — ( £—\1 

- a [f y 


0.015900517 


QJ = 


, f + /f+ 

J gap / J 


l. 000402074 


Q\ = 


-fc / .pc 
J gap / J 


1.07868/8 


Qi = 


' f gap / / 

(fc/f+y-n C +/C c 


1.581883299 


Q2 = 


2.8355305 


94 = 


-ct/[{cmn 2 \ 

■ -ctB 2 /(c+y 


14.697323(20) 


i?4~ = 


7.336774(10) 


r& = 




3.67866(5) 


rs = 


9&/91 


26.041(11) 


no = 


910 /gj 


284.5(2.4) 
4.44(6) x 10 3 


ri2 = 


n 4 = 


gu/gt 


8.43(3) x 10 5 


V 3 = 


-c^B/(c 2 y 

-C 4 -5 2 /(^-) 3 + 3^ 


33.011(6) 


V 4 = 


48.6(1.2) 



Table 8: Critical exponents and universal amplitude ratios for the two-dimensional Ising 
universality class. See App. [3| for the definitions. 
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The constant B^° can be expressed in terms of invariant amplitude ratios: 

= R x (Rt) {1 ~ S)/2 = 0.592357(6) x 10" 5 , (4.10) 

where we have used the numerical results of Table |8|. Moreover, by using Eq. (|4.5| ), we 
obtain 

Bf = c log = 4Q + g 4 = 9.349087(13). (4.11) 

To reach the coexistence curve, corresponding to t < and h — 0, one should perform an 
analytic continuation in the complex t-plane |||7| . The spontaneous magnetization is related 
to the complex zero zq = \zo\e~ t7Tl3 of B(z) |7|], where 

\ Zo \ 2 = R+ = -CtB 2 /(Ctf = 7.336774(10). (4.12) 

Therefore, the description of the coexistence curve is related to the behaviour of B(z) in the 
neighbourhood of zq. 

The functions B(z) and f(x) give equivalent representations of the equation of state. 
Indeed, they are simply related by 

z - 5 B(z) = B^f(x), z = \zo\x~P. (4.13) 



4.2 Parametric representations 

In order to obtain a representation of the critical equation of state that is valid in the whole 
critical region, we need to extend analytically the expansion ( f4.8p to the low-temperature 
region t < 0. For this purpose, one may use parametric representations, which implement 
in a simple way all scaling and analytic properties |8HiO|. One may parametrize M and t in 
terms of R and 9 according to 

M = moR^e, (4.14) 

t = R(i-e 2 ), 

h = h RP 5 h(6), 

where h and m are normalization constants. The variable R is nonnegative and measures 
the distance from the critical point in the (t, h) plane; the critical behaviour is obtained 
for R — > 0. The variable 9 parametrizes the displacements along the lines of constant R. 
The line 9 = corresponds to the HT phase t > and h = 0, the line 9 = 1 to the 
critical isotherm t = 0, and 9 = 9 , where 9q is the smallest positive zero of h(9) to the 
coexistence curve T < T c and h — *■ 0. Of course, one should have 9 > 1. The regularity 
properties of the critical equation of state require h{9) to be analytic in the physical domain 
< 9 < 9 . This is at variance with what happens for the scaling functions f(x) and B(z), 
that are nonanalytic for x — ► oo and z — > oo respectively. This fact is important from a 
practical point of view. Indeed, in order to obtain approximate expressions of the equation 
of state, one can approximate h(9) with analytic functions. The structure of the parametric 
representation automatically ensures the correct analytic properties of the equation of state. 
Note that the mapping (|4.14 ) is invertible only in the region 9 < 9i where 



2 = — ^ = -• (4.15) 
' 1-2/5 3 v ; 
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Thus, the physically relevant interval < 9 < 9q must be contained in the region 9 < 9i, 
and thus we should have 9 < Q\. In practice, since 9i is a singular point of the mapping, 
it is important that 9\ — 9 is not too small. As we shall see, all our approximations satisfy 
this condition. 

The function h(9) is odd in 9, and is normalized so that h(9) = 9 + 0(9 3 ). Since 
M = C^ht^"' for M — > 0, t > 0, this condition implies m = C + h . Following Ref. J7J, we 
then replace h by a single normalization constant p in such a way that we can write 

z = p 9 (l-9 2 )~\ (4.16) 

B(z(9))=p(l-9 2 Y" 5 h(9). (4.17) 

In the exact parametric equation the value of p may be chosen arbitrarily: clearly the physical 
function B(z) does not depend on it. However, if we adopt an approximation for h(9), as 
we will do, the dependence of B(z) on p is not eliminated. One may then choose p to obtain 
an optimal approximation. 

From the function h(9) one may calculate the scaling functions f(x), using the relations 

(4.18) 




Ml) 

and all universal amplitude ratios involving zero-momentum quantities, such as the n-point 



susceptibilities (see, e.g., Ref. O] for a list of formulae) 



4.3 Polynomial approximations for h{&) 



In order to construct approximate parametric representations, we consider a systematic 
approximation scheme based on polynomial approximations of h{9) 0, i.e. 



fe-i 



h^\p,9) = 9 + Y, h ^+i(pW 



2i+l 



(4.19) 



i=l 



The coefficients h 2 i+i(p) are functions of p, 7, /3, and are obtained by matching the small-z 
expansion of B(z) to 0(z 2k ^ 1 ), cf. Eq. (|4.17 ). This kind of approximation turned out to be 
effective for the determination of the critical equation of state of three-dimensional Ising-like 
systems [0,|TTJ, and was generalized to models with Goldstone singularities [37]. In order 
to optimize p for a given truncation h^ k \p,9), we use a procedure based on the physical 
requirement of minimal dependence on p of the resulting universal function 



Bto{p,z)=p(\-PY»hV>{p,0). 



(4.20) 



One may indeed prove |TT] that for any truncation k there exists a solution p^ independent 
of z that satisfies a global stationarity condition, i.e. 



dB^(p, z) 
dp 



0. 



(4.21) 



P=Pk 
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hW 




p 


2.065 


2.027 


2.018(3) 


e? -el 


0.183 


0.177 


0.173 


Tp, 


*3. 67866(5) 


*3. 67866(5) 


*3. 67866(5) 


r 8 


24.413 


♦26.041(11) 


♦26.041(11) 




249.11 


277.1(2) 


♦284.5(2.4) 


r 12 


3513.7(2) 


4066(4) 


4215(49) 


B{\z \/2) 


1.9621 


1.9666 


1.9670 


B(\z \) 


37.160 


41.655 


42.808 


B(2\z \) 


431786 


538946 


569182 




0.4225xl0~ 5 


0.5279(8)xl0~ 5 


0.557(10) xl0~ 5 


Ri 


7.879 


7.558(2) 


7.47(3) 


R x 


7.967 


7.434(4) 


7.23(7) 


u 2 


48.565 


44.41(3) 


42.7(6) 


V3 


28.009 


28.756(6) 


29.2(2) 



Table 9: Universal amplitude ratios obtained by taking different approximations of the 
parametric function h{9). The reported "errors" are only related to the uncertainty of 
the corresponding input parameters. Numbers marked with an asterisk are inputs, not 
predictions. 



for all values of z. 

As input parameters we use the coefficients of the small- z expansion of B(z), i.e. the 
estimates of coefficients r 2 j obtained by TM+CFT and reported in Table [7[ 

In Table |9| we report the universal amplitude ratios derived from truncations correspond- 
ing to k — 3, 4, 5. We will not report the results for k = 6 obtained using the available 
estimate of r\ 2 , because the relatively large uncertainty on r%2 induces a very large error in 
the results of the k = 6 truncation. We only mention that results with k = 6 are perfectly 
consistent with those obtained from the k = 5 truncation. This can be inferred from the 
fact that the estimate of r 12 obtained using h^(9), r 12 — 4215, is very close to the central 
value of the TM+CFT estimate, i.e. r 12 = 4200(740). As we shall see, we will obtain a much 
better estimate of r\ 2 in Sec. fO. 



The results of Table g are not stable as k increases, showing a systematic drift up to 
k = 6, where the large uncertainty does not allow a meaningful comparison. We observe 
that the results for the universal amplitude ratios , R~l , R x , U 2 and v% effectively converge 
towards their precise estimates reported in Table |^. It is also reassuring that the difference 
between the exact value and the estimate obtained using is of the order of the variation 
of the estimates with changing k. 



4.4 Improved approximations from constrained polynomials 

Although the polynomial approximations we presented in the previous Section are substan- 
tially consistent with the known results for the amplitude ratios, they do not provide an 
accurate approximation of the equation of state. The convergence appears rather slow, 
probably requiring the knowledge of a considerably larger number of coefficients r 2 j to sub- 
stantially improve the results. Here, we will present an improved approximation scheme that 
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k 


Pk 


el 




2 


2.01116 


1.15278 


1 


- 0.2084086* 


3 


2.00770 


1.15940 


1 


- 0.2156756> 2 - 0.03 9 4 036* 4 


4 


2.00770 


1.16441 


1 


- 0.2193886* 2 - 0.04179K9 4 - 0.0134886* 6 


5 


2.00881 


1.16951 


1 


- 0.2223896* 2 - 0.0435476* 4 - 0.0148096* 6 - 0.0071686* 8 



Table 10: Polynomial approximations of h(6) using the global stationarity condition for 
various values of the truncation parameter k, cf. Eq. ( 14.251 ). The reported expressions 
correspond to the central values of the input parameters. 



significantly increases the precision of the results. 

The approximation scheme can be improved by constructing constrained polynomial 
approximations of h(8) that take into account the large- z asymptotic behaviour of B(z): 



B(Z) = fin 



(4.22) 



where the value of is reported in Eq. (|4.10| ). We consider constrained polynomial ap- 
proximations of the form 

k— 1 

h^(p, 6) = 6 + Y: h 2t+l (p)6 2l+l + h 2k+1 (p)9 2k+1 , (4.23) 
i=i 

where the coefficients h2i+i{p) with i < k are determined as before, by matching the small- z 
expansion of B(z) to 0(z 2k ~ l ), while h 2 k+i(p) is fixed by requiring that 

= p 1 - 5 h { - k \p, 1) = 0.592357(6) x 10~ 5 . (4.24) 

It follows 

k— 1 

W(P) = P 5 ^B^ - 1 - E h 2 i + i(p). (4.25) 



In this approximation scheme the free parameter p can be still determined by requiring 
the global stationarity condition ( ^4.21| ). This nontrivial property is essentially due to the 



fact that the constraint fl4.24|) is linear in the coefficients /i2i+i- It can be proved by using 



arguments similar to those employed in App. C of Ref. |TT] to show the global stationarity 
condition ( [4.2 1| ) for the approximation scheme ( |4.19| ). In Table ID[ for k = 2,3, 4, 5, we report 
the polynomials (9) obtained by using the global stationarity condition to fix p, and the 
central values of the input parameters Fq°, r$, rg, rio- Note the stability of the coefficients 
of the polynomials with k and that the size of the coefficients decreases with the order of the 
polynomial. The results for some universal quantities are reported in Table |11]. They are in 
much better agreement with the exact results than those obtained without constraint. 

In Fig. [I] we show the scaling function B(z) obtained from TS k \p,6) for k = 2,3,4,5. 
The convergence is good, indeed their differences are not visible in Fig. [I]. This allows us to 
determine B(z) for all real z > with a relative precision of at least a few per mille (the least 
precision is found around z ~ \z§\ ~ 2.71). This fact is not trivial since the small- z expansion 
has a finite convergence radius which is expected to be |z | = (-R4 )^ 2 — 2.71. Therefore, the 



20 











h<® 


p 


2.011 


2.008 


2.008 


2.009(1) 


I 


0.181 


0.174 


0.169 


0.164(2) 


r 6 


3.929 


*3.67866(5) 


*3.67866(5) 


*3.67866(5) 


r 8 


27.585 


26.932 


*26.041(11) 


*26.041(11) 


r 10 


297.25 


292.89 


292.89 


*284.5(2.4) 


r 12 


4425.2 


4385.6 


4385.6 


4443(16) 


fU 


84387 


84029 


84029 


84305(79) 


B(\z \/2) 

V 1 u \l } 


1.9798 


1.9690 


1.9675 


1.9672 


B(\zo\) 

VI u IV 


44.930 


44.335 


44.146(2) 


44.05(3) 


B(3|z |/2) 

V I u 1 / / 


8442.2 


8432.7 


8429.4 


8427.7(5) 


B(2|z |) 

V 1 U|/ 


604619(6) 


604548(6) 


604524(6) 


604511(7) 


5(3|z |) 

V 1 « 1/ 


2.63497(3)xl0 8 


2.63496(3) xlO 8 


2.63495(3) x 10 s 


2.63495(3) x 10 s 




*0.592357(6)xl0- 5 


*0.592357(6)xl0- 5 


*0.592357(6)xl0~ 5 


*0.592357(6)xl0~ 5 




0.021893 


0.021375 


0.021198(3) 


0.02110(3) 


S 2 °° 


9.3987 


9.2611 


9.2611 


9.286(7) 


i2 4 


7.458 


7.396 


7.371 


7.355(5) 


i? x 


7.602 


7.172 


7.002(2) 


6.90(3) 


1/2 


45.918 


42.358 


40.76(2) 


39.6(3) 


"3 


28.328 


29.201 


29.837(9) 


30.5(2) 



Table 11: Universal amplitude ratios obtained from the constrained polynomial approxima- 
tions ( \j.23j ) of the parametric function h(8). The reported "errors" are only related to the 



uncertainty of the corresponding input parameters (they are reported only if they are larger 
than the last figure shown). Numbers marked with an asterisk are inputs, not predictions. 

determination of B(z) on the whole positive real axis from its small- z expansion requires an 
analytic continuation, which is effectively performed by the approximate parametric repre- 
sentations we have considered. We recall that the large- z limit corresponds to the critical 
isotherm t = 0, so that positive real values of z describe the HT phase up to t — 0. Note 
also the good agreement of the results for B%° (see Table |TT]), i.e. the next-next-to- leading 
coefficient of the large- z expansion of B(z), with the precise estimate given in Eq. ( |4.11|) . 
The convergence of the polynomial representations at the coexistence curve is slower. 



This can be seen by looking at the estimates reported in Table |TT| for universal amplitude 
ratios involving quantities related to the coexistence curve, such as Rj[ , R x , and U 2 , and 
comparing them with the corresponding known results reported in Tables ^. They appear to 
(monotonically) converge toward the correct results. The rate of convergence worsens when 
quantities with more and more derivatives with respect to h are involved in the amplitude 
ratio, as it can be already seen by comparing the results for R£ and L^- 

Fig. ^ shows the scaling function f(x) as obtained by the truncations k = 2,3, 4, 5. The 
accuracy of the determination of f(x) can be inferred from the convergence of the curves 
with increasing k and the comparison with the known behavior for x — > —1 and x — > +oo. 



Indeed, for x — >• — 1 we have 



f(x) = b f (l + x) + 0[(l + x) 2 \, (4.26) 



0.69511778. 

it v 
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Figure 1: The scaling function B(z) versus z. We report the curves obtained from the 
constrained polynomial approximations ( 4.2c ) with k = 2,3,4,5. Their differences are not 
visible. 



and for x 



-oo 



f(x)=fS X* + 0(xT-* i ), 

/°° = R- 1 = 0.14752994.... 



(4.27) 



This shows that f(x) is determined with a precision of a few per cent in the whole region. 
This happens also in the large-x region, which corresponds to the HT phase, and therefore to 
z <C 1 in F(z), essentially because f(x) is normalized at the coexistence curve, i.e. x — — 1, 
where our approximation is worse. For x > 0, the error on f(x) increases from to 2%, the 
relative error on R x . 

Using the results of Table |TT| we also obtain 



r 12 = 4.44(6) x 10 3 , 
n 4 = 8.43(3) x 10 5 , 
B™ = 0.0211(2). 



(4.28) 
(4.29) 
(4.30) 



Note that the above-reported estimate of r 12 is perfectly consistent with the result obtained 
by TM+CFT, i.e. r 12 = 4.20(74) x 10 3 , but much more precise. 
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Figure 2: The scaling function f(x) versus x. We report the curves obtained from the 
constrained polynomial approximations ( 4.23Q with k = 2,3, 4, 5. 



5 Conclusions 

Let us briefly summarize the main results of this paper, which presents a determination of the 
equation of state of the two-dimensional Ising model in the whole (t, h) plane. The starting 
point of the analysis was the determination of the 2rz-point couplings constants in the HT 
phase which are related to the expansion of the free energy in powers of the magnetization. 
These quantities were then used as input parameters for a systematic approximation scheme, 
which allowed us to obtain an accurate determination of the equation of state in the critical 
regime in the whole (t, h) plane. 

The determination of the coefficients r2j of the small-magnetization expansion of the 
Helmholtz free energy is in itself an interesting and nontrivial problem. We addressed this 
problem by using TM techniques. In order to overcome the typical problem of all TM 
calculations, i.e. the fact that they are constrained to rather small values of the lattice size 
in the transverse direction (in our case we could reach L\ = 24 as maximum value) we used 
a two-step strategy: 

1] In the first step we studied the model at very small values of j3 so as to have small values 
for £ (whose maximum value was £ = 3.350288.., which was reached for the largest 
value of (3 that we studied: (3 = 0.37) and thus large values of the ratio We 
then extrapolated to the thermodynamic limit the TM results by using a new effective 
iterative algorithm. This algorithm is rather interesting in itself and could be of general 
utility for people working with TM methods. 



23 



2] Second, we used a combination of standard RG results and a set of high-precision nu- 
merical informations on the zero-field magnetic susceptibility to construct the scaling 
functions for the derivatives of the free energy involved in the construction of the r 2 j 
coefficients. With these scaling functions we could obtain the critical-limit values of the 
r2j coefficients with rather small errors even if the input data for this continuum- limit 
extrapolation correspond to rather small values of /3. 

The expansion of the free energy in powers of the magnetization in the HT phase was then 
used as a starting point to construct approximate representations of the equation of state 
that are valid in the critical regime of the whole (t, h) plane. We considered a systematic 
approximation scheme based on polynomial parametric representations, devised to match the 
known terms of the small-magnetization expansion and the large-magnetization behaviour 
of the Helmholtz free energy. A global stationarity condition was used in order to optimize 
the polynomial approximation. This approximation scheme can be improved systematically 
by considering higher and higher order polynomials. It is only limited by the number of 
known terms in the small-magnetization expansion of the free energy. The knowledge of this 
expansion up to 10th order allowed us to obtain an accurate determination of the critical 
equation of state. We indeed obtained the scaling function B(z), cf. Eq. ( |4.7|) , for all real 
z > with a relative precision of at least a few per mille, and the scaling function f(x), 
cf. Eq. fl4.1|) , with a precision of a few per cent in the whole physical region x > — 1. The 
approximation scheme is systematic, thus this precision can be improved by a more accurate 
knowledge of the small-magnetization expansion of the free energy in the HT phase. 

The method that we used to reconstruct the equation of state from the small-magnetization 
expansion of the Helmholtz free energy is general and can be applied to other statistical mod- 
els. We mention that similar methods have been successfully applied to the three-dimensional 
Ising ffTliPH an d XY universality f37|,[38| classes, leading to accurate determinations of the 



critical equation of state and of the universal ratios of amplitudes that can be extracted from 
it. 
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A Notations for the amplitudes 



Universal amplitude ratios characterize the critical behaviour of thermodynamic quantities 
that do not depend on the normalizations of the external (e.g. magnetic) field, order param- 
eter (e.g. magnetization) and temperature. Amplitude ratios of zero-momentum quantities 
can be derived from the critical equation of state. Beside the amplitudes of the n-point 
functions Xnif) defined in Eq. (|3.1|) and the corresponding ones C~ in the low-temperature 
phase, we consider amplitudes derived from the singular behaviour of the specific heat 

C H = A ± \n(l/\t\) ) (A.l) 

the spontaneous magnetization on the coexistence curve 

M = B{-tf. (A.2) 

We complete our list of amplitudes by considering the second-moment correlation length 

k^rtr, (a.3) 

and the true (on-shell) correlation length, describing the large-distance behaviour of the 
two-point function, 

u. = or- ( A -4) 

One can also define amplitudes along the critical isotherm, e.g., 

x = or 7//35 , (a.5) 

£ = f a \h\- v W, (A.6) 
£ gap = flJhr^. (A.7) 



B Finite-size scaling of the free energy 



Before starting the discussion on the FSS in the Ising model, let us stress that the analysis 
reported in this appendix is only a straightforward application of results which are already 
well-known in the literature. They can be found for instance in ]39| and are based on the 



results of [ifj and on the exact solution of the Ising model on a finite lattice obtained by 
|T2"|] . We decided to write it down in this appendix all the same so as to make this 

to simplify the 



Kaufman 



39 



paper as self-contained as possible. We shall use the same notations as 
comparison. This means in particular that we shall use the letter R to denote the finite size 
of the lattice in the transverse direction, which is denoted in the rest of the paper with L\. 

The aim of this appendix is to obtain the explicit form of the FSS of the free energy 
of the Ising model on a rectangular lattice of size R x oo, which is exactly the geometry 
we considered in our TM work. To obtain this result we shall work in the framework of 
the S-matrix approach to 2d integrable models and shall in particular use the so-called 
Thermodynamic Bethe Ansatz (TBA). In the TBA one looks at the theory defined on an 
infinitely long cylinder with circumference R. The only free parameter of the theory is 
r = mR where m is the lowest mass of the model. The goal is to extract the behaviour of 
the free energy E (R) and of the lowest mass m(R) (and possibly, in models more complicated 
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than the Ising one, of higher massive states) as a function of R. The geometry of the TBA 
is exactly the one which we have in our TM setting. The parameter r is the ratio L/£ in our 
notations and the only thing that we must require to compare our findings with the TBA 
analysis is that we should be in the "scaling region," i.e. that £ ^> 1 so that we may neglect 
the lattice artifacts. 

The TBA results depend on the entries of the S-matrix. They are in general very compli- 
cated. Thus it is usually impossible to find the exact FSS functions for any value of r (one 
has usually to resort to perturbative expansions for small or large values of r) . However, the 
thermal perturbation of the Ising model is so simple (the S matrix is simply S = — 1 as a 
consequence of the fact that the model can be mapped to a free massive field theory) that 
in this case the explicit expression for any value of r can be obtained. 



B.l Free energy at h = 

The TBA prediction for the FSS behaviour of the free energy, in the HT phase of the 2d 
Ising model at h — 0, isQ 

F{R) = F(oc) - ^ (B.l) 



where, following the standard TBA machinery (see p. 668 of Ref. |39|]), c(r) is given by 

6 r°° 



TV 2 JO 



/ d9 rcosh( 9) ln( 1 + e- rcoshe ) . (B.2) 
Jo 



The integral can be exactly evaluated and gives: 

c(r) = ^- 2 J2 [ -^—K 1 (kr), (B.3) 
71 k=i K 

where K\ is the modified Bessel function of the second kind. 

The normalization of Eq. ( |B.1|) could seem strange, but it is chosen such that in the 
r — > limit (i.e. at the critical point) the function c(r) exactly becomes the central charge 
of the Ising model c = 1/2. The idea behind this choice is the so called c-theorem of 
Zamolodchikov which states that the FSS of any 2d integrable model can be parametrized 
by Eq. ( |B.1| ), where c(r) is a "running central charge," which in general interpolates between 
two critical points and, in this case, between the central charge of the Ising model and the 
value c(oo) = which is the appropriate one for a massive theory. 

The limit in which we are interested is r 3> 1. In this limit, we can use the asymptotic 
expansion for the Bessel functions and Eqs. ( B.l ), ( B.3 ) become: 



F w= F <»>-S£ e -""^< 1+ -> (B ' 4) 



where the dots stand for the o(l) terms. 



3 Notice that in |59[ the authors study, instead of the free energy F(R), the quantity Eq(R) that is related 
to F{R) by E (R) = R{F(R) - F(oo)}. 
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B.2 Derivatives of the free energy at h = 



Since we are interested in the derivatives with respect to h of the free energy, we must try 
to extend Eq. ( |B.4| ) in the h ^ plane. It is quite reasonable to assume that in ( |B.4| ) the 
h dependence can only be hidden in the mass term. Moreover, we know from Refs. fill, 42 1 
that this dependence (for small values of h) is well described by an analytic even function of 
h, i.e. 

m(h) =m(l + a h 2 + b h A + ...), (B.5) 

with a and b unknown constants. 

Inserting this expression in Eq. (|B.4|) and performing the derivatives we find 



X2n(i?)-X2n(oo) = R 2 ^' 2 (g 2n>1 ( 1 / R)e~ mR + g 2n , 2 (l/i?)e" 2miJ + g 2n , 3 ( l/R)e~ 3mR + . . .) (B.6) 

where the g 2n ^(l/R) may be expanded in positive powers of 1/R (whose coefficients could 
be in principle computed as functions of a, b, and m). This is the expression quoted in the 
text as Eq. (fTLj). 



B.3 A test of the iterative algorithm for the infinite- volume ex- 
trapolation 

In order to further check the IA ( 3.14 ) for the infinite-size extrapolation, we apply it to a 
test function of the type (P-6|). We consider the function 



E(£,R) = 1 + 



7r <9 4 c(r) 



6R 2 <9/i 4 

where c(r) is the function defined in Eq. (P-3| ) and 



(B.7) 



h=0 



R 
1 



(l + h 2 + K 



(B.8) 



Clearly, for R — > oo, E(£, oo) = 1. We compute E(£, R) for finite values of R in the typical 
range of our TM calculations, i.e. 3 ^ R/£ ^ 7, and apply the IA (|3.14j) to determine 
the E(£, oo). The comparison with the exact value gives an idea of the effectiveness of the 
method. Table [T| shows the results of the IA fl3.14|) for £ = 4, analogously to Table |3] for 
_F 4 . We report the results for various iteration levels (up to four); the level zero corresponds 
to the original data. In order to obtain from Table [12] an estimate E(£, oo), we follow the 
same strategy used in Sec. |3.1.2|, i.e. we consider the largest value of R and four iterations 



of the IA. The error is estimated from the residual dependence on R of the results obtained 
with four iterations of the algorithm. We obtain E(£ = 4, oo) = 1.000000(1), which is in 
perfect agreement with the exact number. Analogous results are obtained for other values 
of £. This confirms the effectiveness of the procedure we used to perform the extrapolation 
to the thermodynamic limit, and the reliability of the uncertainty we considered. 
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R 








£(3) 




12 


1.002919604 
















13 


1 


0038256116 
















14 


1 


0041103404 


1 


0042408301 












15 


1 


0040272017 


1 


0040459912 












16 


1 


0037413786 


I 


0041444428 


1 


0041 1 1 3946 








17 


1 


0033581849 


1 


00486621 18 


1 


004030441 








18 


1 


0029428048 


I 


0083034836 


1 


0039525998 


1 


0020057776 




19 


1 


0025338931 


{) 


9766831 739 


1 


0052032241 


1 


0040258799 




20 


1 


0021528510 


{) 


9969431 584 


o 


989031 3495 


1 


0040423723 


1.0040425080 


21 


1 


0018100907 





9987411366 





9989162382 





9951663087 


1.0040259105 


22 


1 


0015092173 





9993480457 





9996572961 





9997173548 


0.9981748024 


23 


1 


0012497886 





9996258634 





9998603960 





9999370740 


0.9999482199 


24 


1 


0010291081 





9997722737 





9999353998 





9999793168 


0.9999893714 


25 


1 


0008433668 





9998559476 





9999675467 





9999916598 


0.9999967552 


26 


1 


0006883533 





9999063496 





9999827009 





9999962157 


0.9999988812 


27 


1 


0005598763 





9999378507 





9999903522 





9999981547 


0.9999995914 


28 


1 


0004540024 





9999580850 





9999944245 





9999990581 


0.9999998462 



Table 12: Results of the IA pJD applied to tie test function E(£ = i,R), cf. Eq. flBwf) . 



C Determination of 

In this appendix we present the determination of e^. We will obtain it by analyzing the 
scaling corrections to the free energy in the presence of a magnetic field h on the critical 
isotherm t = 0. The scaling corrections for h — > have been determined in ||. Using the 
notation of 0, we have, cf. Eq. (84) of ref. [g: 

F(0,h)= f b +A) |/i|T§(l + 4 ib |/i|ti+4 )1 |/i|ii+4 )2 |/ i |i+4 )3 |/ i |i + •••). (C.l) 
Each of these amplitudes has a precise physical meaning. Let us look at them in detail: 

fb denotes the bulk contribution to the free energy and can be obtained from the exact 
solution of the Ising model on the lattice at the critical point. Explicitly /& = -G + 
| log 2, where G is Catalan's constant. 

A l f is the amplitude of the singular part of the free energy. It can be evaluated exactly in the 
framework of the S-matrix approach to the model. Its value is A l j = 0.9927995... @,0. 

Af b is the first correction (proportional to h 2 ) due to the bulk part of the free energy. 



It can be related to the constant term Dq = —0.104133245. . . [|3] appearing in the 
expansion in powers of t of the susceptibility at h — 0. Its value is A l j b = D /(2A l f) = 
-0.0524442... 

A l f i is due to the TT, T 2 and T 2 irrelevant operators in the Hamiltonian. This amplitude 
turns out to be compatible with zero || . 

A l f 2 is due to the b t h 2 term in u t . 
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A l f 3 is due to the eth 2 term in Uh- 

In principle one could use the TM data to estimate all these constants. In practice 
this procedure works only for the first unknown amplitude. All the higher ones are then 
"shadowed" by the first. 

The problem with is that it appears at a rather high level. In Q we were able to fix 
exactly the amplitudes only up to A l f V The first unknown one was A l j 2 which was then 
estimated numerically, obtaining 

0.020 < A l f 2 < 0.022 . (C.2) 

It was impossible to give any reliable estimate for A l f 3 . The main progress of the present 
paper with respect to that analysis is that, thanks to the exact calculation of bt performed 
in [|S],[n]] we are now in the position to estimate exactly also A l f 2 - A direct calculation gives 

4 2 = 4^, (C3) 

where A l j was define above, Eq is given by Eq. ( |3.21| ), and A l E is defined by the singular 
behaviour of the internal energy, i.e. E int (t = 0, h) = \dF/d(3 = E hnlk + A l E h 8 / 15 + ... 
Numerically, A l E = 0.58051... H,g§. Substituting in Eq. (Q, we find 

A f , 2 = 0.0210115... (C.4) 

in perfect agreement with the estimate of ^Jj. 

Substituting in Eq. ( |C.1| ), we may now estimate numerically the amplitude A l ^ 3 , which 
is related to by: 

Af,s = ^e h . (C.5) 

A standard application of the fitting procedure discussed in J§], using as input data those 
reported in Table 10 of || gives 

e h = -0.00727(15). (C.6) 
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